Interfacial and wetting properties of a binary point Yukawa fluid 
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We investigate the interfacial phase behavior of a binary fluid mixture composed of repulsive point 
Yukawa particles. Using a simple approximation for the Helmholtz free energy functional, which 
yields the random phase approximation (RPA) for the pair direct correlation functions, we calculate 
the equilibrium fluid density profiles of the two species of particles adsorbed at a planar wall. 
We show that for a particular choice (repulsive exponential) of the wall potentials and the fluid 
pair-potential parameters, the Euler-Lagrange equations for the equilibrium fluid density profiles 
may be transformed into a single ordinary differential equation and the profiles obtained by a 
simple quadrature. For certain other choices of the fluid pair-potential parameters fluid-fluid phase 
separation of the bulk fluid is observed. We find that when such a mixture is exposed to a planar 
hard-wall, the fluid exhibits complete wetting on the species 2 poor side of the binodal, i.e. we 
observe a thick film of fluid rich in species 2 adsorbed at the hard-wall. The thickness of the wetting 
film grows logarithmically with the concentration difference between the fluid state-point and the 
binodal and is proportional to the bulk correlation length of the intruding (wetting) fluid phase. 
However, for state points on the binodal that are further from the critical point, we find there is no 
thick wetting film. We determine the accompanying line of first-order (pre-wetting) surface phase 
transitions which separate a thin and thick adsorbed film. We show that for some other choices of 
repulsive wall potentials the pre-wetting line is still present, but its location and extent in the phase 
diagram is strongly dependent on the wall-fluid interaction parameters. 

PACS numbers: 68.08.Bc, 61.46.-w, 05.20.Jj, 05.70.Np 



I. INTRODUCTION 

When charged nano-particles (colloidal macroions) are 
dispersed in a neutralizing medium the electrostatic po- 
tential between the particles is screened by the counteri- 
ons and is much shorter-ranged than the bare Coulomb 
interaction between the two identical particles [l], 
Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, for 
example, predicts that the screened electrostatic contri- 
bution to the effective interaction potential between the 
particles takes the form Q: 



Q* A exp (-Ar) 
47re r r 



(1) 



where Q* is the re-normalized charge, e r is the dielec- 
tric constant of the solvent and A -1 is the Debye screen- 
ing length which determines the thickness of the double 
layer of opposite charge surrounding each colloidal par- 
ticle. The DLVO theory also includes a hard-core re- 
pulsion, that takes into account the size of the colloids, 
and dispersion (van der Waals) attraction. In many sit- 
uations the latter plays a minor role and the effective 
particle-particle interaction can be modeled by a hard- 
core, repulsive Yukawa potential. If one considers the 
limit of high charge and/or a low density of particles, 
one finds that the particles infrequently come into con- 
tact and the effect of the hard-core on the properties of 
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the fluid is small. Under these circumstances one may 
argue [3| that Eq. |T]), applied for all separations r, is a 
good zeroth-order model to describe a charged colloidal 
suspension. Indeed, Hynninen and Dijkstra ,3| showed 
that the phase diagrams of one-component hard-core re- 
pulsive Yukawa particles could be mapped to those of a 
point Yukawa system for sufficiently large charge on the 
particles. 

The Yukawa potential does, of course, arise in a variety 
of physical situations. For example, the point Yukawa 
pair potential may also be used to model the interactions 
between micron sized dust particles in a charge neutral 
plasma - so called dusty plasmas Q. The binary point 
Yukawa model has also been used in a recent simulation 
study of glassy dynamics in low temperature-low density 
(classical Wigner) glasses @- 

In a previous study [|| , we investigated the bulk struc- 
ture and phase behavior of the two component point 
Yukawa fluid, i.e. a binary mixture in which all the po- 
tentials between the particles are repulsive Yukawas. We 
found that for certain choices of the parameters in the 
pair potentials, and at sufficiently high densities, the 
mixture separates into two fluid phases, one of which is 
composed predominantly of particles of species 1 and the 
other phase is predominantly of species 2. In this paper 
we investigate the inhomogeneous mixture adsorbed at a 
single planar wall. Using a simple approximate density 
functional theory (DFT), we calculate the equilibrium 
fluid density profiles and thermodynamic quantities rel- 
evant for investigating surface phase behavior; the latter 
are surprisingly rich given the simplicity of the model and 
DFT. 
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In Ref. [6| we showed that by making comparisons 
with results from the accurate hypernetted-chain (HNC) 
approximation the simple random phase approximation 
(RPA) for the pair direct correlation functions Cy 



a, 



with i,j — 1,2, in the Ornstein-Zernike equation 
yields a reasonably accurate approximation for the bulk 
fluid radial distribution functions Qijir) at large and in- 
termediate r, provided that the fluid is at a state-point 
fairly close to the fluid-fluid binodal (specifically, within 
the Lifshitz line [6|). However, the RPA was found to 
be inaccurate for determining the fluid correlation func- 
tions far from the coexistence region, where one must 
use the more reliable HNC or other related bulk theories 
@j H 5 E3] • We also found that the RPA provides a good 
approximation for determining thermodynamic quanti- 
ties such as the bulk fluid Helmholtz free energy. Indeed, 
the phase diagram resulting from the RPA is in fairly 
good agreement with that obtained from the HNC Q. 

In the present work we use a simple approximate DFT 
that, in bulk, generates the RPA for Cij{r) to investigate 
the inhomogeneous mixture. Based on the experience 
from studying the bulk fluid mixture Q we expect that 
the DFT should predict with reasonable accuracy the 
fluid density profiles and adsorption behavior for state 
points that are close to the binodal. In such cases, where 
a fluid close to fluid-fluid phase coexistence is adsorbed 
at a planar wall, it is possible for there to be wetting 
of the wall by the coexisting phase, for some choices of 
the wall-fluid interaction parameters [ll|. Specifically if 
a mixture rich in particles of species 1 is at a state point 
close to the coexistence line, then we may observe wet- 
ting of the wall by a thick film of the coexisting phase 
that is rich in species 2. This is indeed what we observe, 
for several different choices of purely repulsive wall po- 
tentials. Furthermore, we find that the present system 
may exhibit a surface phase transition from a thick to a 
thin adsorbed wetting film and we examine the depen- 
dence of this transition on the specific form chosen for 
the wall potentials. 

Some of the inspiration for the present study comes 
from the particularl y el egant DFT studies of wetting pi- 
oneered by Sullivan [12j for a one component fluid, where 
the fluid pair potential featured a hard-core plus an at- 
tractive Yukawa tail. The DFT that Sullivan used was 
similar to the one we use here. The hard-core repulsion 
(absent in the present system) was treated by making 
a local density approximation, and the Yukawa attrac- 
tion treated in a simple mean-field approximation. In 
a subsequent study of the corresponding binary mixture 
adsorbed at a hard-wall with exponential wall-fluid at- 
tractive potentials, Telo da Gama and Evans [HI showed 
that for a particular set of wall-potential and fluid pair- 
potential parameters [3] , the Euler-Lagrange (EL) equa- 
tions arising from the minimization of the density func- 
tional yield a single ordinary differential equation (ODE) 
that can be integrated directly to determine the equilib- 
rium fluid density profiles. The benefit of this approach is 
that not only can the density profiles and surface tension 



be calculated easily without the need for sophisticated 
numerical schemes, but also that the criteria for differ- 
ent types of film adsorption and locating wetting transi- 
tions could be established directly [H, [H[ • We apply an 
equivalent method to the two component point Yukawa 
model, treated within the present DFT, and show, in 
Appendix A, that for a particular set of fluid-fluid and 
wall-fluid interaction parameters, we may also derive a 
single ODE that determines the equilibrium density pro- 
files. However, in the present case the restriction on the 
parameters precludes fluid-fluid phase separation so that 
wetting must be investigated by solving the coupled EL 
equations. 

This article proceeds as follows: In Sec. Ill Al we in- 
troduce the binary point Yukawa fluid model and recall 
some key results from Ref. [6|. Sec. IIIBI describes the 
inhomogeneous fluid situation and defines the external 
(wall) potentials under consideration. In Sec. Ill CI we 
briefly describe the approximate DFT that generates the 
RPA and outline the general DFT approach to calculat- 
ing adsorption behavior. We present numerical results 
for the fluid density profiles, adsorption, and for the lo- 
cation of the pre-wetting transition line, for various wall 
potentials in Sec. IIIII Finally in Sec. IIVI we discuss our 
results and draw some conclusions. 



II. BACKGROUND 
A. The Model Fluid and its Bulk Properties 

We consider a binary fluid composed of particles where 
the pair potentials between particles of species i and j are 
given by: 



4>ij (r) 



Mije exp(— Ar) 



4/T 



Ar 



(2) 



where e denotes the overall energy scale, My > are the 
(dimcnsionless) species specific interaction magnitudes, 
and A is an inverse decay-length. Following our previous 
work [|| , we allow the inter-species interaction magnitude 
A/12 = (1 + 5)t/MixM22 to depend on a non-ideality 
parameter d. For point ions immersed in a medium with 
inverse screening length A, 6 = 0. However, as discussed 
in Q , one can contemplate situations in charged colloidal 
mixtures where charge renormalization leads to effective 
potentials with 5 =/= 0. In the examples that follow, we fix 
the dimensionless temperature T* = ksT/e = (/3e) _1 = 
1, where ks is Boltzmann's constant, and we set M\\ = 1, 
M22 = 4 and S = or 0.1. In Ref. [6j], we showed that 
for an ideal mixture, 6 = 0, the fluid does not phase 
separate, while for any S > there is phase separation 
at some, sufficiently high, densities. Note that in Eq. @ 
there is an additional factor 1/Att that was not used in 
the definition of the pair potentials in Ref. Q . 

The bulk phase diagram for the model fluid with 
S = 0.1 calculated using the RPA is displayed in Fig. [T] 
Details of how this phase diagram is obtained are given 
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FIG. 1: The bulk phase diagram for the binary Yukawa fluid, 
for the set of pair-potential parameters Mn — 1, M22 = 4 
and S = 0.1, calculated using the RPA. p b is the total density 
and £2 is the concentration of species 2. At sufficiently high 
densities the fluid demixes. The two-phase region is bounded 
by the binodal (solid line), which meets the spinodal (dotted 
line) at the critical point (o). The straight tie-lines connect 
coexisting state-points with pressures (3PX~ 3 = 150 to 500, in 
increments of 50, and then from 500 to 1900 in increments of 
100 (from bottom to top). The arrows mark three paths which 
terminate at the binodal (on the species 2 poor side), which 
are at fixed total densities A: p b X~ a = 30, B: p b X~ i = 24, 
and C: p b X~ i — 18. The inset shows a magnification of these 
paths. In subsequent figures we display results for the fluid 
density profiles at a hard-wall for various X2 along these paths. 



in Ref. [6|. This mixture exhibits a two-phase fluid-fluid 
coexistence region which is bounded by the binodal (solid 
line). Within the binodal is the spinodal (dotted line), 
which is the line in the phase diagram at which the RPA 
predicts that the compressibility and the bulk correlation 
length diverge. The binodal and the spinodal meet at a 
single critical point (circle) in the phase diagram. Within 
the RPA, scaling the interaction potentials <fiij(r) by a 
factor A results in a scaling of the total density, p b , by a 
factor A -1 . Thus, in Fig. Q]the coexisting densities are 
larger by a factor of 47r than those in the phase diagram 
displayed in Fig. 4 of Ref. f(|. Note that the latter figure 
compares RPA results for the binodal and spinodal with 
those obtained from the HNC. 

The general theory for the asymptotic decay, r — > 00, 
of the bulk fluid total pair correlation functions hij (r) — 
9ij ( r ) — 1) states that for a system of particles interacting 
via short ranged potentials, such as in the present system, 
ha (r) may be obtained from the following expression jy, 

m 



rhij(r) 



exp(ig„r), 



(3) 



where the summation is made over contributions from 
the poles in the upper half of the complex plane of the 
function D(q), which is a non- linear combination of the 



Fourier transforms of the pair direct correlation functions 
Cij(<l) [ILIIII- The asymptotic decay, r — ► 00, of hij(r) is 
determined by the pole q n with the smallest imaginary 
part an and A% is the amplitude. 

In Ref. ,6] it was shown that the asymptotic decay, 
r — > 00, of hij(r) for the present system, obtained using 
the RPA, is determined by one of two purely imaginary 
poles. The location of these poles for the rescaled pair 
potentials, Eq. |2]), is given by: 



On 



2AT* 



M ± 



^M% + Ms)+\ 2 , (4) 



where p b is the total density, Mn = (1 — g)Mn + X2M22 
and M 5 = 4(1 - x 2 )x 2 M 11 M22{2 + 6)6 [ff . For the case 
when 6 = 0, it can be shown that there exists only one 
pole with ccq > A. For the case 6 > 0, there is a second 
pole, «q < A. Of these two poles, the one with the 
smaller an determines the ultimate asymptotic decay of 
hij(r) and also defines the bulk fluid correlation length, 



that £ > A -1 . 



r hff A (r) ~ A- ex p(~ a o r), v — > 00. Note 



B. Wall-Fluid Potentials 

We investigate adsorption of the model fluid at a pla- 
nar substrate. A single planar wall exerts the potential 
Vi(z) on particles of species i, where z is the Cartesian 
axis perpendicular to the wall. We consider here three 
different wall-fluid potentials. These are all based on the 
hard-wall, which is infinite for z < and zero for z > 0: 



Vi(z) = 



00 z < 







z > 0, 



(5) 



for i = 1,2. Note that this potential is not dependent 
on any energy or length scale parameters. The hard-wall 
potential may be modified by adding an extra interaction 
term for z > 0. This contribution may be repulsive, 
attractive, or a combination of both. Here we restrict 
consideration to purely repulsive potentials. Adding an 
exponential tail maintains the infinite step discontinuity 
at the origin: 



Vi(z) 



Ait exp(— Xz) 



z<0 
z > 0, 



(6) 



where the amplitudes A, > for both species of parti- 
cles. This choice of external potential corresponds to the 
case where the hard-wall carries a (screened) charge that 
is uniformly distributed over the surface. Furthermore 
this potential allows us to recast the EL equations as an 
ODE [12]. We also consider a Yukawa wall potential, 
similar to the inter-particle interaction potential, of the 
form 



VAz) = \ °° Z < 

1 Ajeexp(— Xz)/Xz z> 



z < 
z > 0, 



(7) 
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where both amplitudes Aj > 0. Note that the Yukawa 
wall potential removes the infinite step discontinuity in 
Vi(z) at the origin, replacing it by a smooth increase. In 
both cases we set the wall potential decay length scale 
equal to the inter-particle interaction length scale A -1 , 
in order to simplify the model. This choice therefore pre- 
cludes any phenomena that may arise from a competition 
between different length scales in the potentials [16[ . We 
require that the density profiles decay to their bulk values 
far from the wall, 



lim pi{z) 



P b i, 



for i = 1,2 and determine the adsorption of species i, 
Ti, which in the present planar geometry is given by the 
integral 



r, 



Az{p l (z) - p\). 



(9) 



In regions of the phase diagram where the density profiles 
vary smoothly as a function of changing state-variable 
(e.g. the concentration or total density), then we also 
expect Ti to vary smoothly. Conversely, where the den- 
sity profiles vary discontinuously as a function of a state- 
variable, discontinuities will occur in both Ti and I^, and 
will signal surface phase transitions. 



C. Implementation of Density Functional Theory 

Here we present a brief description of the DFT a 
proach. For a more complete account see e.g. Refs 
E3; (3 • For a fluid composed of v different species of par- 
ticles the thermodynamic grand potential is a functional 
of the set of one-body density profiles {pi(r)}, i = 1 . . . v; 



V[{pi}} = F[{ Pi }} - J dr ft (r)( Mi 



Vf*(r)), (10) 



where ^[{pi}] is the intrinsic Hclmholtz free energy func- 
tional, [ii is the chemical potential of species i and V^ ext (r) 
is the external potential acting on species i. Minimizing 
the grand potential functional with respect to variations 
in the density profiles, one obtains a set of v EL equa- 
tions that may be solved simultaneously to obtain the 
equilibrium fluid density profiles: 



SF[{ Pi }} 



+ Vr t (r). 



(11) 



The intrinsic Helmholtz free energy functional may be 
separated into a sum of two contributions: ^[{pi}] = 
-Fid[{Pi}] + -Pex[{Pi}]- The first is the ideal-gas term and 
the second is the excess contribution due to the particle 
interactions. The ideal gas term is: 



Fid[{ Pi }] =Y,l 3 ~ 1 f dr ft (r)(ln(Afo(r)) 



1), (12) 



where Aj is the (irrelevant) thermal de Broglie wave- 
length of particles from species i. In the present study 
we consider a binary mixture [y = 2) , and employ a sim- 
ple mean-field approximation for the excess part of the 
Helmholtz free energy: 



2 

E 



J J drdr'piWp^rO^-dr-r'l). 



This functional generates the RPA approximation for 

,-(|r - r'l) = 



(8) the pair direct correlation functions: 



-0 



-/3<fe(|r-r'|) [191,12c 



' S P i(r)S Pj (r') ~ 

Since we only consider external potentials that vary in 
one Cartesian direction z, Eqs. (JSHZJ), the fluid density 
profiles also only vary in the z-direction. Using Eq. (|13p , 
we can write Eq. (fTTJ) as follows 



Vi{z) = (iid,i(z) + 



E 

3=1 



dr'^Gr-r'l)^'), (14) 



where Pid,i{z) — fi~~ x hi(A| Pi{z)) is the ideal-gas contri- 
bution to the chemical potential of species i. Given the 
external potentials and the boundary conditions in Eq. 
(|8j, we may solve the EL equations (fl4|) to obtain the 
equilibrium fluid density profiles. Generally, the solu- 
tions must be obtained numerically, for example, by using 
a simple Picard iterative scheme. However, for a partic- 
ular choice of (exponential) wall and particle interaction 
parameters, one can show that the EL equations may be 
transformed to yield a single ODE that can be integrated 
once to give an explicit quadrature for the fluid density 
profiles. The derivation is outlined in Appendix A. 



III. RESULTS OF CALCULATIONS 

We have performed calculations for a variety of bulk 
state points and for various choices of the wall-fluid and 
inter-particle potential parameters. Representative re- 
sults are presented here. 

A. Exponential wall and 5 — 0.0 



For the particular case when S — 0, and the exponen- 
tial wall potential © with A2 = y^j^Ax, we can use 

the approach described in Appendix [A] to calculate the 
fluid density profiles. Recall that for 6 — the fluid 
does not demix at any density [6(. From the analysis 
of Appendix A we find that the density profiles, Pi{z) 
are both monotonic functions of z, and the densities at 
the wall are determined by Eq. (|A11[) . For the given 
set of pair potential parameters, the wall- fluid potential 
magnitude, A\, determines the behavior of the fluid at 
the wall. For A\ 0, the wall-fluid potential is equiva- 
lent to a hard-wall and we find strong adsorption of both 
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FIG. 2: The density profiles p%{z) of a binary point Yukawa 
fluid at an exponential wall, Eq. ([6]), calculated by solving 
the ODE, i.e. from Eq. ((SJ and (|A3|) . The solid line is 
pi{z) and the dashed line is p2{z). The parameters for the 
fluid pair potentials are Mn = 1, M22 = 4, and <5 = 0, and 
the repu lsive expo nential wall parameters are A\ = 1 and 
A 2 = i/M 22 /Mii^i = 2. The bulk fluid at z —> 00 has 
p b A _3 = 30 and X2 = 0.5. The density profiles exhibit strong 
adsorption at the wall but decay rapidly and monotonically 
to the bulk values. 



species at the wall. Recall that for a hard wall, the sum 
of the contact densities is equal to the bulk fluid pres- 
sure divided by ksT and the pressure is large for these 
mixtures. For A\ ~ 1 there is again strong adsorption 
at the wall, but the value of the contact density is de- 
creased. For A\ ^> 1 then pi(0) oc e~^ tAl resulting in a 
depletion of both species at the wall and monotonically 
increasing profiles. For all state points each density pro- 
file returns to its bulk value over a distance ~ £, the bulk 
correlation length. In Fig. [2] we display a typical set of 
density profiles showing such behavior for the case when 
the total bulk density p b = p\ + p\ = 30A 3 and the bulk 
concentration X2 = p\j ' p b = 0.5. 

In order to understand the adsorption at the wall it is 
important to consider the contributions to the free energy 
that arise from the fluid-fluid and fluid-wall interactions. 
We first consider the simplest case, A\ — *■ 0. The in- 
troduction of the hard-wall results in the total absence 
of any fluid behind the wall, i.e. pi(z < 0) = 0, for 

1 = 1,2. Therefore, any particles close to the wall have 
reduced their fluid-fluid interaction energy, which in turn 
favors strong adsorption of both species at the wall. Since 
particles of species 2 have a larger (repulsive) interac- 
tion energy, the decrease in the free energy is greater for 
species 2, which results in stronger adsorption of species 

2 at the wall. However, as A± increases there is an addi- 
tional free energy contribution for particles close to the 
wall due to the increased wall-fluid repulsive interaction. 
If A\ is increased sufficiently then this energy contribu- 
tion negates the benefit from particles adsorbing at the 
wall and the adsorption is reduced. As A\ is increased 



FIG. 3: The density profiles pi(z) for the point Yukawa mix- 
ture with 8 = 0.1 at a hard- wall, Eq. ([5]), calculated for the 
fluid at state points along path A in Fig. [1] The density 
profiles of species 1 are shown in the main figure and the 
species 2 profiles are shown in the inset. For all concentra- 
tions X2, up to and including bulk coexistence, both density 
profiles show strong adsorption at the wall, but decay to the 
bulk densities over a distance ~ £, the bulk fluid correlation 
length. Results are shown for X2 = 0.0011, X2 = 0.0022 and 
X2 = 0.0033 ~ X2, C ocx, the coexistence concentration for this 
total density, p 6 A~ 3 = 30. The density profiles for the fluid 
at coexistence are those with the dashed lines. 



still further then ultimately both species become depleted 
at the wall, since the wall-fluid interaction becomes the 
dominant contribution to the free energy. This rather 
simple behavior is not unexpected. Recall that the bi- 
nary fluid does not demix for 5 = 0, so we do not observe 
any wetting or layering behavior. 



B. Hard wall and 8 = 0.1 

For mixtures with 8 — 0.1 the coupled integral equa- 
tions (|14p must be solved numerically, rather than us- 
ing the approach described in Appendix A. In this sub- 
section we describe the behavior of the fluid at the hard- 
wall, with potentials given by Eq. ([5]). We find that the 
adsorbed fluid exhibits different types of behavior, de- 
pending on the state-point. However, one common fea- 
ture is that there is always strong adsorption at the wall. 
The contact densities, Pi(0), for both species of parti- 
cles can be rather large, and this effect is particularly 
marked in the case of the particles of species 2, which 
are always favored by the hard-wall, and stems from 
the high pressures in these systems. For state-points 
far away from bulk coexistence we find that for both 
species the density profiles exhibit a thin adsorbed layer 
at the wall, and decay to their bulk values over a distance 
~ £ = 1/o.q > A -1 . For all state-points on the right side 
of the binodal in Fig.[TJ i.e. rich in species 2, this type of 
decay persists up to and including the binodal. However, 
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FIG. 4: As in Fig. [3j except here the density profiles are 
calculated along path C in Fig. [1] at constant total density 
p b X~ A = 18. For concentrations close to the binodal the den- 
sity profiles exhibit a thick wetting film rich in species 2 in- 
truding between the bulk fluid and the wall. The thickness 
of this film diverges logarithmically where 
a;2,cocx — 0.04207 is the bulk coexistence value. The density 
profiles displayed are for concentrations from X2 = 0.022 to 
X2 = 0.042 ~ a;2,cocx, in increments of 0.004. The dashed-line 
density profiles are those for 22 — a;2.cocx where the wall is 
almost completely wet by the phase rich is species 2. 



to the left of the binodal in Fig. [T] i.e. for state-points 
poor in species 2, we find that in some instances the 
adsorbed layer grows much thicker on approaching the 
binodal, while for some other state points only a thin 
adsorbed layer remains. In order to illustrate this, we 
display results for density profiles calculated along three 
different paths, increasing xi towards the binodal, at con- 
stant total densities. These three paths A, B, and C, at 
successively lower densities, are displayed in Fig. [TJ 

On path A, with total density p b \~ 3 = 30, wc find that 
for all concentrations, the density profiles decay to their 
bulk values over a finite distance ~ £, so that the adsorp- 
tions, Ti and T2, remain finite up to and including the 
coexistence state point. A series of density profiles along 
path A are displayed in Fig. [3l Note the strong adsorp- 
tion and high contact densities at the wall. Along path 
C, at constant density p b \~ 3 = 18, we find that the den- 
sity profiles are radically different for state-points close 
to coexistence. Instead of there being a thin (finite) ad- 
sorbed layer at the wall, we find a thick film of fluid rich 
in species 2 adsorbed at the wall. As X2 — > £"2,cocx, where 
^2,coox is the concentration at the binodal, the thickness 
of the wetting film increases and ultimately diverges, i.e. 
the adsorption T2 — * +00 as X2 —> £2,cocx- This behavior 
is termed complete wetting. A number of density pro- 
files along path C are displayed in Fig. [H indicating the 
growth of the wetting film as X2 — ► a^.cocx- Note that the 
density profiles of species 1, Pi(z), are non-monotonic. 
As X2 increases, the contact density pi(0) reduces, and 
a minimum develops in the profile at Xz < 0.5. The 



growth of the wetting film is accompanied by increasing 
depletion of species 1, and Ti — > —00 as X2 —> a^.cocx- 

Recall that for a one-component fluid, exhibiting 
liquid-gas phase coexistence, one observes that for a bulk 
gas with density near to coexistence, a thick wetting film 
of the liquid may be adsorbed at the planar substrate 
or wall. In the case when all the potentials are short 
ranged, it can be shown that the thickness of this film, 
I, diverges as I ~ — Iq In \p g — p b \ where p g is the den- 
sity of the gas at coexistence. The amplitude Iq depends 
on the relative ranges of the wall-fluid and the fluid-fluid 
potentials [ll|, [HI- Provided I is large, the adsorption 
at the wall T is proportional to I: T ~ (pi — p g )l, where 
pi is the density of the coexisting liquid. In an entirely 
analogous way, in the present two component system, as 
X2 approaches coexistence at constant p b , the adsorption 
of species 2, T 2 , is given by 



r 2 



p 2 ' a ) In \x 2 - x 2 



(15) 



where p h 2 is the bulk coexisting density of the (wetting) 
phase (3 rich in species 2 and p^'" is the same quan- 
tity in phase a, poor in species 2 (l6| . For the present 
wall potentials, Eqs. I©-©, we expect the length scale 
^0 = £,w, the bulk correlation length of the fluid phase 
rich in species 2 that wets the wall. The adsorption calcu- 
lated from the DFT via Eq. ([9]) is plotted for state-points 
along path C in Fig. [5^ where we also compare with the 
asymptotic result in Eq. (fT5|) with Iq = £ w = 1 /a^ , where 
Q!q is the imaginary part of the pole that determines the 
asymptotic decay of hij (r) in the wetting phase, obtained 
from Eq. ©. Note that since cvq depends on the total 
bulk density, it is different in the two coexisting phases. 
The asymptotic formula for Ti, equivalent to Eq. (fl~5j) . 
provides an equally good fit to the DFT results for this 
quantity. 

It is important at this stage to understand why we get 
wetting for 6 — 0.1, when for S — we did not. For S = 
0.1 there is an added energy cost for particles of species 
1 and 2 to mix, which in bulk drives phase separation 
at sufficiently high densities. It is the combination of 
this effect and the tendency of the fluid to adsorb at the 
wall that induces wetting. The high densities of both 
species at the wall promotes local phase separation; in 
this case creating a region rich in species 2 and poor in 
species 1 which increases as the state-point approaches 
the binodal. 

We now consider the fluid interfacial behavior along 
path B, which has constant total density p b A~ 3 = 24.0. 
As we increase X2, we find initially that the density pro- 
files decay rapidly to the bulk values over a short distance 
~ £, similar to the density profiles on path A. However, 
as we increase the concentration X2 further, we find a 
discontinuous change in the density profiles which leads 
to a discontinuous change in the adsorptions IV Increas- 
ing X2 towards the coexistence value, we find that the 
profiles are similar to those on path C and that T2 di- 
verges in a manner similar to that described for path C 
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FIG. 5: As in Fig. [3j except here the density profiles are 
calculated along path B in Fig. [T] at constant total density 
p b \- 3 = 24.0. The concentrations are x 2 = 0.0091, 0.0095, 
0.0096, 0.0097 and 0.0098. We observe a discontinuous change 
in the density profiles as X2 is changed continuously. The thin 
to thick adsorbed film transition (pre- wetting) occurs between 
X2 = 0.0095 and 0.0096. The coexistence concentration at this 
total density is a;2,cocx = 0.0113. 



as X2 — > a;2,cocx- This discontinuous change in the ad- 
sorption denotes a point on the pre- wetting phase transi- 
tion line, which is a line in the phase diagram that sepa- 
rates regions with thick and thin adsorbed films [2l[. In 
Fig. [5] we display a number of density profiles calculated 
on path B, showing the discontinuous change on cross- 
ing the pre-wetting line. In Fig. [6)d we display the DFT 
results for T2, calculated from the density profiles, and 
compare with the asymptotic result in Eq. (|15p . where 
once again io = Cw = 1/ &o * s ^ ne bulk correlation length 
of the wetting phase, calculated from Eq. ([4J. 

In order to establish the location of the pre-wetting 
transition line, it is necessary to calculate the grand po- 
tential, f2. In the vicinity of the pre-wetting transition we 
find that there are two branches of solutions which min- 
imize the grand potential - one branch corresponding to 
a thin adsorbed film and the other to a thick adsorbed 
film. We compute the value of SI for each set of pro- 
files to establish which set is the global minimum of the 
grand potential. In practice, we calculate f2[{pi}], scan- 
ning along lines of constant total density in the phase- 
diagram, for both increasing and decreasing concentra- 
tion, x%. By plotting fl versus x-2 we determine the con- 
centration where these two branches intersect (where the 
grand potentials of the thick and thin film branches are 
equal) which is the pre-wetting concentration at the par- 
ticular total density. Since the gradient of the minimum 
grand potential changes discontinuously and the adsorp- 
tions, Ti, jump discontinuously, this is a first order phase 
transition. By repeating this procedure for different total 
densities, we are able to map out the location of the pre- 
wetting line in the phase diagram. We find that it extends 
tangentially from the binodal, as discussed in Ref. [22l |. 



u 1 



80.0 



70.0 



60.0 



50.0 



40.0 



30.0 




-10.0 



W-Voex"^ 1 



FIG. 6: The adsorption of species 2, F2, (points) at a hard- 
wall plotted versus In \x2 — a;2,cocx|, the logarithm of the dif- 
ference between the species 2 concentration and the value at 
bulk coexistence, calculated a) along path C in Fig. [T] which 
is at total density p b \~ 3 — 18.0, and b) along path B, where 
p b X~ i = 24.0. The solid line is the asymptotic result in Eq. 
(115[) . where lo = £ w is the bulk correlation length in the phase 
wetting the wall, calculated from the RPA - see text. In b) 
we observe a jump in T2 at \n\x2 — £2,cocx| — —6.3, which 
corresponds to the pre-wetting transition. 



and terminates in a pre-wetting critical point, where the 
difference in adsorption Ti between the thin and thick 
adsorbed film vanishes. The hard- wall pre-wetting line 
(long dashed line) is displayed in Fig. \7\ As the inset 
shows, the line lies rather close to the bulk binodal and 
extends over a narrow range of densities. 



C. Repulsive Exponential and Yukawa walls 

Since the hard-wall potentials do not have any ad- 
justable parameters, the location of the pre-wetting line 
is fixed for a given set of fluid-fluid interaction parame- 
ters. However, if we consider a hard- wall augmented by 
repulsive exponential or Yukawa tails, Eqs. © and ([7]), 
this introduces two independent parameters (the ampli- 
tudes Ai) which may be varied. Studying various com- 
binations of these parameters, we find that the occur- 
rence of wetting, and the existence and location of a pre- 
wetting line is strongly dependent on the magnitudes of 
A\ and A^. For A\ = A2 — > we recover the hard- 
wall behavior. For combinations of parameters where 
Ai ~ \JMa with i = 1,2, we find that with both the 
exponential and Yukawa walls, there is a pre-wetting line 
that is not drastically removed in the phase diagram from 
the location of the hard-wall pre-wetting line - see Fig. [7] 
for examples with A\ = 1 and A2 = 2. If we increase 
the ratio A2/A1, we find that the pre-wetting line moves 
down the binodal, towards the bulk critical point. As it 
approaches the critical point, there are indications that 
the wetting transition may change from a first order to 
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x 2 

FIG. 7: The pre-wetting lines (PWL) for the binary point 
Yukawa fluid at a hard-wall and at hard-walls augmented by 
repulsive exponential and Yukawa tails. The solid line is a 
portion of the bulk binodal from Fig. [T] The main figure 
is a magnification of the inset, which displays the spinodal, 
the bulk critical point (circle) and the hard-wall PWL. All the 
PWLs descend tangentially from the binodal and terminate in 
a critical point, denoted by the symbol •. The hard-wall PWL 
meets the binodal at p b \~ z = 28.2 and the critical point is 
located at the density p^~ s = 22.5 and concentration X2, c = 
0.012. The PWL for the hard-wall with repulsive exponential 
tail Eq. ((6} is for the case with parameters A\ — 1.0, Ai — 2.0. 
Much further down the binodal is the PWL for the hard-wall 
with repulsive Yukawa tail, Eq. (0, with parameters Ai = 
1.0, A2 = 2.0. The PWLs for the exponential and Yukawa 
walls are much shorter than the hard-wall PWL, although by 
varying A\ and A2 the location and extent of the PWL may 
be varied - see text. 

a continuous (critical) wetting transition, although we 
have not investigated this issue in detail - when locat- 
ing the wetting transition numerically, it can be difficult 
to discriminate between a first order wetting transition 
with a very short pre-wetting line and a true continuous 
wetting transition. If A2/A1 is increased even further 
the wetting transition and all wetting behavior may be 
pushed onto the opposite branch of the binodal, i.e. any 
wetting is by the phase rich in species 1. Conversely, if 
the ratio A2/A1 is decreased we find that the pre-wetting 
line moves up the binodal creating a larger complete wet- 
ting regime. In all cases where there is complete wetting, 
we confirmed that the increase of T2 as the binodal is ap- 
proached is given by Eq. (fT5|) with Iq = irrespective 
of the magnitudes of A\ and A 2 . 



IV. CONCLUDING REMARKS 

Using a simple mean field DFT, we have investigated 
the interfacial behavior of a two component point Yukawa 
fluid adsorbed at a planar wall. For a restricted set of 
the parameters in the model, the EL equations reduce to 
a single ODE that can be integrated to yield a simple 



expression, Eq. (|A9[) . for the density profiles. However, 
the restriction on the parameters required to make this 
simplification limits the method to cases where the fluid 
does not exhibit fluid-fluid phase separation, i.e. 6 = 0, 
and therefore precludes the study of wetting behavior. 
This scenario differs from that in the binary mixture of 
hard-core plus attractive Yukawa tail particles adsorbed 
at a hard-wall augmented by exponential attractive tails 
where the equivalent mixing rules lead to different classes 
of wetting behavior [l3[ . In the general unrestricted case, 
we must obtain the density profiles numerically by solv- 
ing either the coupled pair of integral EL equations in 
Eq. (fT4|) . or for the exponential wall potential Eq. ([6]), 
the coupled pair of differential equations in Eq. (|A2|) . 
For the model fluid mixture that exhibits phase sepa- 
ration (where the parameter 8 > 0), we find generally 
that particles of species 2 are more strongly adsorbed at 
the wall than are particles of species 1. Of course, one 
could choose the wall potential parameters, by making 
A 2 » Ai in Eqs. © and ([7]), so that the effective at- 
traction between the wall and the particles of species 1 is 
sufficiently strong that this situation is reversed. In the 
more common case, where there is a stronger effective 
attraction between the wall and the particles of species 
2, we find that for some state points on the species 1 rich 
side of the bulk binodal, that a thick wetting film rich in 
species 2 forms at the wall. In particular this occurs for a 
hard wall. We find that in the complete wetting regime, 
the thickness of this wetting film increases logarithmi- 
cally as the bulk concentration approaches its value at 
coexistence. 

Generally we find a first-order wetting transition with 
a pre-wetting line extending out of bulk coexistence for 
the various wall potentials that we have considered. How- 
ever, the location of this pre-wetting line is very sensitive 
to the precise form of the wall potentials and to the val- 
ues of the parameters in the various potentials. It is 
important to note that for all the wall potentials that 
we considered, the decay length of the potentials is A , 
which is the length scale in all the fluid pair-potentials - 
see Eq. @. We have not investigated the influence on the 
wetting behavior arising from varying this length scale. 
We expect such a modification to have a strong influence, 
not only on the existence and location of the pre-wetting 
transition line, but also on the 'rate' of the growth of the 
thick wetting film [c.f. Eq. (|15p ]. In a study of a different 
binary fluid, with short ranged Gaussian potentials, ex- 
hibiting bulk fluid phase behavior that is similar to the 
present system [l6l |20|. it was found that if the wall po- 
tentials are of the form in Eqs. ^ or (O, but with with 
a modified decay length 7 -1 (i.e. with the length A -1 re- 
placed by the decay length 7 -1 in these potentials), then 
the amplitude Iq in Eq. (|15[) is no longer necessarily the 
bulk correlation length of the wetting phase, £ w . In the 
case of the exponential wall-potential Eq. ([6]) , the longer 
of 7 _1 and £ w determines the prefactor l in Eq. (fT5|> 
[TTI . [l|| . In the case of the Yukawa wall-potential J7]) , if 
£iu > 7 _1 , then Iq = £ w . However, when £ w < 7 -1 , then 
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the prefactor 1$ in Eq. (fT5|) is neither £ w nor 7 _1 
We expect the same scenario for the present model fluid. 

One important question to address is how robust are 
the present results. Are the phenomena that we observe 
simply an artifact of using the simple (RPA) DFT? In 
to order address this question we have performed further 
calculations, that we do not describe in detail here, where 
we have studied the interfacial phase behavior, deter- 
mined the location of the pre-wetting line and examined 
the growth of the wetting film using a different (more so- 
phisticated) approximation for the excess Helmholtz free 
energy functional. The functional that we utilized is the 
following: 

F ex [{ Pi }} = F^M}] +5>r / dr^(r) 

E / dr / dr'^(r)^(r')c y (|r-r'|), (16) 

which is obtained by making a Taylor expansion of the 
excess Helmholtz free energy functional around that of 
the uniform fluid with densities {p\} and truncating the 
expansion at second order in Spi(r) = pi(r) — p\ [l|[. 
-Fea;[{Pi}] i s the excess Helmholtz free energy of the uni- 
form (bulk) system, p,? x is the excess chemical potential 
of species i in the uniform system and cy (r) are the pair 
direct correlation functions in the bulk reference fluid far 
from the wall. For all of these bulk fluid quantities, that 
are required as inputs to the theory, we use results ob- 
tained from the HNC theory. We find that for the differ- 
ent wall potentials that we considered above within the 
(RPA) DFT the location of the pre-wetting line obtained 
from the functional [)16p and the RPA functional (fT5|) are 
located very close together, especially when the moder- 
ately small differences in the location of the binodals ob- 
tained within the two theories are taken into account Q . 
Thus we are confident that all of the phenomena that we 
have observed within the simple RPA theory are at least 
qualitatively correct. At present we do not know of any 
computer simulation results for the binary point Yukawa 
fluid that would confirm this assertion. 

Although the point Yukawa model of a fluid is very 
simplistic, and may be inappropriate for modeling the 
wide range of real systems where hard-core effects domi- 
nate the physics, its usefulness lies in its ability to incor- 
porate realistic and relatively complex behavior, without 
the need for an elaborate density functional theory. We 
have shown that with a very simple model and functional 
it is possible to describe rich wetting behavior, which is 
also reproduced in the results of a more sophisticated 
functional. In order to study effects arising from an ex- 
plicit hard-core one would require another (hard-sphere) 
contribution to the functional e.g. 18], and a significantly 
greater computational effort to calculate the equilibrium 
density profiles. This becomes especially important when 
we move away from simple planar or spherical geome- 
tries to situations where the densities vary in two or even 
three dimensions. For example, we have used the binary 



point Yukawa model to study the effective interaction be- 
tween a large colloidal particle and a thick wetting film 
adsorbed at a planar wall, and a fluid-fluid interface [23| ■ 
The use of the simple RPA functional greatly decreases 
the computational cost of calculating individual equilib- 
rium profiles, and of establishing the location and nature 
of phase transitions. 

Interfacial phase behavior similar to that presented 
here, i.e. the existence of a first-order pre-wetting tran- 
sition and logarithmic growth of thick wetting films has 
also been observed in other binary systems of soft-core 
particles namely the Gaussian core and star-polymer so- 
lutions [H, HII- These similarities lead us to conclude 
that the behavior found for the present model should be 
quite generic to binary mixtures of purely repulsive par- 
ticles that exhibit fluid-fluid demixing. 

Finally, we note that since the present DFT treatment 
is purely mean-field, the influence of capillary wave fluc- 
tuations in the wetting film interfaces is neglected. In the 
complete wetting regime fluctuations merely change the 
amplitude of the logarithmic growth from £ to £(l+a>/2), 
where u> is the standard dimensionless parameter mea- 
suring the strength of the fluctuations [ll], [lj| [25| . We 
believe that a treatment taking fluctuations into account 
would arrive at the same prediction of a first order pre- 
wetting phase transition. Where fluctuations are likely 
to be more important is in the tri-critical regime where 
the cross-over from first-order to critical wetting occurs. 
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APPENDIX A: A SINGLE ORDER-PARAMETER 
TREATMENT OF THE MIXTURE. 

We follow closely the derivation in Refs. [H, [H| f° r 
a related model fluid. We begin by setting the external 
potentials to the exponential type, defined in Eq. ([6]). 
Substituting Eqs. ^ and © into the EL equation JTj 
gives 



(Md,i( x ) = Pi - Aieexp(-x) 



(Al) 



f-f 2A 3 



da;' exp(— \x — x'\) P j(x'), 



where x = Xz is a dimensionless length. Note that since 
the density profile is zero for z < 0, as a consequence 
of the hard external potential, the lower integral limit 
becomes 0. Taking two derivatives with respect to x on 
both sides of Eq. (|A1[) we obtain the following expression: 



d 2 pid,t( x ) 
dx 2 



p\dA x ) — im. + x 



2 



Mijepjix). (A2) 
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which is equivalent to that derived in [13] for a mix- 
ture with hard-core pair potentials and Yukawa tails ad- 
sorbed at a hard-wall with exponential tails; for the point 
Yukawa model the hard sphere chemical potential Mh,i is 
replaced by the ideal gas chemical potential, Mid.i- 

Eq. (|A2|) constitutes a pair of coupled equations (i = 
1.2). By substituting one into the other and rearranging, 
it can be shown that these two equations for the fluid 
density profiles may be made independent of each other 
if and only if M 12 = ^M 11 M 2 2 [14]. This corresponds 
to the ideal geometric mixing rule, where 6 = 0. It is 
also necessary for the external potential parameters to 
be related by the rule: A 2 = ^/M^jM^Ai [II]. Using 
these parameters the density profiles for the two different 
species are then related by a simple scaling factor: 



M, 



Mid,2(a0 - M2 = \l (Mid,i(aO - Ml)- 



(A3) 



Having reduced the problem to determining a single func- 
tion, it is possible to rearrange Eq. (|A2|) to obtain a single 
ODE for the single order parameter p\d,i{x) that deter- 
mines the density profiles of both species: Using the local 
Gibbs-Duhcm relation 



pi(x) 



i = 1,2 



(A4) 



where Pid{x) = kBT[pi(x) + P2{x)] is the local ideal-gas 
contribution to the pressure and setting i = 1, we may 
now write Eq. (|A2|) as 



-^Mid.i(z) = Wd,i(x) - Mi + MiieA-^ dpid( f v (A5) 
dx z d/x id4 (x) 

where the total derivative is [c.f. Eq. (|A3|) ]: 



where p\ is the bulk density of species i. By integrating 
both sides of Eq. (|A5|) . and using the boundary condi- 
tions to ensure that the pressure and the chemical po- 
tentials tend to their bulk values as x — > oo, we find that 



d/iid,i(x) 



-i 2 



d.r 



= [pidAx) - Mi] 2 + 2M n eX 3 [pid(x) - p] 
= ^(Mid.i)- (A8) 



This equation provides an implicit relation for the chem- 
ical potential, fJ.id,i( x ), and therefore the density profile 
Pi(x). The final step is to integrate: 



d^id, 



,i(0) ±VW*id~ 



(A9) 



Note that the order parameter p,id,i(x) which follows 
is necessarily a monotonic function of x and the choice 
of sign in (IA9[) depends on whether the order parame- 
ter is an increasing or decreasing function of x. Since 
Pid,i(x) = fi^ 1 \n(K\pi{x)) it follows that the density 
profile p\{x) is also monotonic. In order to calculate an 
explicit solution, we must use a numerical method such 
as the Runge-Kutta method [26] . The only input that is 
required to this equation is pid,i{x = 0), the value of the 
ideal gas chemical potential at the wall. By differentiat- 
ing Eq. (|A1|) . evaluating at x = 0, and substituting back, 
one finds 



dMidq(O) 
da; 



Mid.i + 2Aie - mi- 



(A10) 



This may be combined with Eq. (|A8|) . evaluated at x = 0, 
to yield the following relation for Mid,i(0): 



d 



d 



22 



d 



(A6) 



dMid,i(^) <9Mid,i(a:) S Mudp^aix)' 
The total pressure, p, of the bulk system can be written 



(3p = Ppid + 



[MM)' + 2M l2 p\p b 2 + M 22 (p b 2 ) 



p5 + /4 + ^£ 



A 



M 



22 P\ 



(A7) 



2Ai[wd,i(0) + Axe- mi] =M u A-> id (0)-p]. (All) 

Note that if — > then we recover the contact-density 
sum- rule for a hard- wall; pi(0) + p2(0) = (3p. Thus, given 
the particle pair interaction parameters M\\ and M 22 , 
and the wall-potential parameter A±, one may calculate 
Pi(x) from Eq. (|A9|) . where the solution of Eq. (|A11|) 
is used as input. One may then calculate p 2 {x) from 
Eq. (|A3[) : p 2 {x) is also monotonic. Results from this 
approach are displayed in Sec. IIIIt \. 
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